A SCUBA/Spitzer Investigation of the Far Infra- Red 
Extragalactic Background 
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ABSTRACT 

We have measured the contribution of submillimeter and mid-infrared sources to the extra- 
galactic background radiation at 70 and 160/zm. Specifically, we have stacked flux in 70 and 
160/zm Spitzer Space Telescope (Spitzer) observations of the Canada-UK Deep Sub-millimeter 
Survey 14h field at the positions of 850/im sources detected by SCUBA and also 8 and 24/zm 
sources detected by Spitzer. We find that per source, the SCUBA galaxies are the strongest and 
the 8/zm sources the weakest contributors to the background flux at both 70 and 160/zm. Our 
estimate of the contribution of the SCUBA sources is higher than previous estimates. However, 
expressed as a total contribution, the full 8/zm source catalogue accounts for twice the total 24/im 
source contribution and ~ 10 times the total SCUBA source contribution. The 8/im sources ac- 
count for the majority of the background radiation at 160/im with a flux of 0.87±0.16 MJy/sr and 
at least a third at 70/im with a flux of 0.103±0.019 MJy/sr. These measurements are consistent 
with current lower limits on the background at 70 and 160/im. Finally, we have investigated the 
70 and 160/im emission from the 8 and 24/xm sources as a function of redshift. We find that 
the average 70/zm flux per 24/zm source and the average 160/im flux per 8 and 24/zm source is 
constant over all redshifts, up to z ~ 4. In contrast, the low-redshift half (z < 1) of the of 8/zm 
sample contributes approximately four times the total 70/zm flux of the high-redshift half. These 
trends can be explained by a single non-evolving SED. 

Subject headings: Cosmology: ; Galaxies 



1. Introduction 

Excluding the microwave background, approx- 
imately half of the entire extragalactic back- 
ground radiation is emitted by dust at far infra- 
red (IR) and sub-millimeter (submm) wave- 



lengths fe.g.jFixsen et al. Hl998tlHauser fc Dwek 
20011 iDole et al. Il2006j) . This cosmic IR back- 
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ground (CIB) radiation peaks at a wavelength of 
~ 200/zm, yet compared to the optical, relatively 
little is known about the sources responsible. 

Surveys conducted by the Sub-millimeter Com- 
mon User Bolometer Array (SCUBA) and the 
Max-Planck Millimeter Bolometer (MAMBO) 
over the last decade have directly resolved up 
to two-thirds of the CIB at 850/zm and 1.1mm 
into discrete, high redshift sources (although this 
fraction is uncertain due to the uncertainty in 
measurements of the CIB at these wavelengths). 
Discovery of this population has been extremely 
important since SCUBA galaxies represent the 
most energetic star-forming systems at an epoch 
when the Universe was at its most active. How- 
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ever, the impact this has had on understanding 
the nature of the CIB is relatively minor since at 
these wavelengths, the CIB has 30 to 40 times less 
power than at the peak. 

A recent study using a large s ample of 73 bright 
(>5m Jy) SCUBA sources by IChapman et al 



( 20051 ) indicated that the population contributes 
a mere ~ 2% of the CIB at the peak, with an 
extrapolation of up to ~ 6% for sources down to 
the fainter limit of lmJy However, this work re- 
lied on an assumed spectral energy distribution 
(SED) for the SCUBA sources constrained only 
by a redshift, the 850/im SCUBA flux and a ra- 
dio flux at 1.4GHz. Furthermore, redshifts were 
obtained from optical spectra having identified 
the optical sources with radio counterparts to the 
SCUBA sources. This introduces two selection 
effects. The first causes SCUBA sources with 
z > 3 to be missed by requiring a radio detection, 
the selection function for radio sources falling off 
rapidly at z ~ 3 due to the K-correction. The 
second causes a paucity of sources around z ~ 1.5 
where no emission lines fall within the observable 
wavelength range of their spectra. 

This motivates the first of two main goals 
of this paper. By stacking the flux in 70 and 
160/im MIPS images at the positions of SCUBA 
sources detected in the Canada-United Kingdom 
Deep Sub-millimeter Survey (CUDSS) 14- hour 



field (|Eales et al. I l2000l : IWebb et al. 1 l2003h . 



directly measure the contribution of the SCUBA 
sources to the CIB at wavelengths in the vicinity 
of the peak. 

In addition to the submm surveys, space-borne 
mid-IR surveys conducted by the Infra-red Astro- 
nomical Satellite (IRAS) and the Infra-red Space 
Observatory (ISO) have resolved significant con- 
tributions to the CIB from the shorter wave- 
length si de of the peak (see, for examp le, the re- 
view bv lLagache. Puget fc Dole I [2005( 1. The in- 
troduction of the Sp itzer Space Telescope (Spitzer; 



Werner et al. 1120041 ) means that such surveys can 



be carried out over much wider areas and to much 
greater depths. 

In particula r, the Multi-band Photometer for 
Spitzer (MIPS; [pdeke et al. I [20041 ) has been used 
for a variety of mid- and far-IR s urveys to resolve 



source s contributing to the CIB. iPapovich et al. 



(|2004l ) showed that approximately 70% of the CIB 
at 24/im can be resolved into IR galaxies with flux 



> 60/xJy. In contrast, iDole et al. I (|2004f ) found 
that at the longest two MIPS wavelengths, 70 
and 160/im, only 20% and 10% of the CIB can 
be directly resolved into distinguishable sources 
brighter than 3.2 and 40mJy respectively. How- 
ever, the error on these fractional quantities is very 
large since the absolute flux of the background at 
160/im is currently unknown to a factor of ~ 2 and 
at 70/im, the uncertainty is even larger. 

A major problem with attempting 
rectly resolve sources in deep 70 and 
MIPS surveys is source confusion due 
large instrument point spread function 
This problem can be circumvented by measur- 
ing the 70 and 160/im MIPS flux at the po- 
sition of objects selected in other wavebands 
for which there are already accurate positions. 
This stacking technique has been successfully 
used by several authors with SC UBA data that 
also suffer from confusion 



to di- 
160/im 
to the 
(PSF). 



Peacock et al. 



20001: ISerieant et al. 1120041: iKnudsen et al 112 005: 



Dve et al. 1 120061: I Wang. Cowie fc Barger 1 120061 ). 
Dole et al. I (|2006h stack MIPS flux at the posi- 
tions of 24/im sources with fluxes > 60/tJy to find 
that they represent the bulk of the CIB at 70 and 
160/im respectively (see Section PO]) . These con- 
tributions are investigated as a function of 24/im 
source flux and, based on external studies of the 
redshift distribution of MIPS 24/im sources, the 
authors conclude that the majority of the radia- 
tion must be emitted at z ~ 1. 

This provides the second main motivation for 
the present paper . We extend the analysis of 
Dole et al. I (|2006l ) in two ways Firstly, we ad- 
ditionally measure the contribution to the CIB 
at 70 and 160/im from 8/im sources observed 
with Spitzer's In fra-red Array Camera (IRAC; 



Fazio et al. 



20041 ). Secondly, we investigate how 
the contribution from the 8 and 24/im populations 
varies with redshift, using photometric redshifts 
established for th ese sources in our earlier work 
(|Dve et al. Il2006h . 



This paper is set out as follows. In the follow- 
ing section we describe the data. Section[3]outlincs 
our stacking procedure. Our results are presented 
in Section 01 followed by a summary and brief dis- 
cussion in Section [5J 
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2. Data 

Coverage of the CUDSS 14h field in this paper 
comprises 850/im SCUBA observations as well as 
data acquired with both Spitzer's IRAC and MIPS 
instruments. The photometric redshifts of the 24 
and 8/im sources used later were determined from 
ground-based U, B, V, I and K photometry as well 
as IRAC 3.6/ mi and 4.5/xm obse rvations. We refer 



the reader to lDve et al. for a full account 



of the determination of these redshifts. 

2.1. SCUBA data 

The SCUBA catalogue contains sources ex- 
tracted from 63 hours worth of 850/im data taken 
on 20 different nights over the period from March 
1998 to May 1999 at the James Clark Maxwell 
Telescope (JCMT). The 850/im map of the ~ 
7' x 6' survey region was composed by combin- 
ing several jiggle maps at different base positions. 
Each jiggle map was observed for approximately 
one hour with a 64-point pattern (to ensure full 
sampling), nodding JCMT's secondary mirror and 
chopping by 30" in right asc ension. We refer the 



reader to lEales et al~T(|2000l ) for more specific de- 
tails of the data reduction. 

The source list used for the sta cking is that 
compi led from the 850/im data by IWebb et al. 



( 20031 ). consisting of 23 sources above a 3tr detec- 
tion threshold within the 41 arcmin 2 SCUBA map. 
The average 3<r sensitivity limit of the sample is 
3.5mJy. 20 of these sources lie within the MIPS 
70/im coverage and 22 within the 160/im coverage. 

2.2. Spitzer Space Telescope Data 

The Spitzer observations discussed in this pa- 
per were obtained as part of the Guaranteed Time 
Observing program number 8 to image the ex- 
tended Groth strip, a 2° x 10' area at a ~ 14 h 19 m , 
5 ~ 52°48' (J2000) with IRAC and MIPS. In the 
present work, we have limited the stacking to a 
small section of the extended Groth strip that fully 
contains the CUDSS 14 hour field. This ensures 
a self-consistent comparison between the SCUBA, 
24/im and 8/im source stacking. Of this section, 
96% of the CUDSS 14 hour field falls inside the 
70/im coverage and 91% inside the 160/im cov- 
erage. The south-east corner of the 70/im data 
and the south-east and north-east corners of the 
160/im data have either poor or no coverage and 



these are masked out in all analyses throughout 
this paper. Figure [T] shows the images. 

The 24 and 8/im source catalo gues used for 



stacki ng were first presented in lAshbv et al. 
( 20061 ) and we refer the reader to this work for 
a detailed account of their creation. Both cata- 
logues cover a slightly larger area than the original 
CUDSS 14 hour field but are contained by the 70 
and 160/im image sections described above. The 
24/im sources cover an area of 49 arcmin 2 and all 
lie above a 5cr point source sensitivity of 70/iJy. 
The 8/im sources cover 59 arcmin 2 and lie above 
a 5cr point source sensitivity of 5.8/iJy. There 
are a total of 177 24/im sources that lie within 
the MIPS coverage at 70/im and 171 within the 
160/im coverage. Of the 8/im sources, 801 and 773 
lie within the MIPS coverage at 70/im and 160/im 
respectively. 

The MIPS 70 and 160/im images were observed 
in scan map mode with the slow scan rate. The 
data were processed with the Spitzer Scie nce Cen- 



tre (SSC) pipeline ([Gordon et al. 1 120051 ) to pro 



duce images with flux measured in MIPS instru- 
mental units. These were converted to units 
of mJy/arcsec 2 using the calibration factors 14.9 
mjy/arcsec 2 per data unit for the 70/im data and 
1.0 mJy/arcsec 2 per data unit for the 160/im data. 
Note that these are 5 — 10% smaller than those 
quoted for the MIPS-Ge pipeline in the MIPS data 
handbook (version 3.2.1, 6 Feb 2006 releas^H) since 
the SSC pipeline is completely independent. The 
pixel size in the 70 and 160/im images is respec- 
tively 9.85" x 9.85" and 16.00" x 16.00". For com- 
parison, the FWHM (full width at half max) of 
the instrumental PSFs are ~ 16" at 70/im and 
~ 40" at 160/im (measured by fitting to the central 
Gaussian component of the PSF) . The data have a 
5ct point source sensitivity of lOmJy at 70/im and 
60m Jy at 160/im. 

The error images output by the current version 
of the SSC pipeline are only an estimate of the 
true error and do not accommodate for the full 
range of effects exhibited by the MIPS detectors. 
Also, the MIPS image data are covariant as a re- 
sult of pixel interpolation and rebinning carried 
out during pipeline construction of the mosaics. 
This covariance must be quantified and incorpo- 
rated into the stacking analysis that follows. For 



1 see \hitp://ssc. spitzer. caltech. edu/mips/dh\ 
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these reasons, we generated our own error data. 

To obtain variance images, we made the as- 
sumption that the error in the flux of a given pixel 
is inversely proportional to the square root of the 
number of times the pixel has been scanned. Of 
course, pixels containing bright sources will have 
an additional contribution from Poisson noise, but 
since our data have very few bright sources and 
since we investigate their effect on the stacking by 
removing them, this is not a concern. The vari- 

14 h 18 m 10 s 18 m 00 s 50 s 40 s 30 s 20 s 




14 h 18 m 10 s 18 m 00 s 50 s 40 s 30 ! 




Fig. 1 — MIPS observations of the CUDSS 14h field 
at 70/im (top) and 160^im (bottom). Sources detected 
with > 3<r significance are circled. Blank regions are 
masked areas of poor or no coverage. 



ance images are then the inverse of the coverage 
maps scaled to have a standard deviation of unity. 
Indeed, we recover a perfect Gaussian distribution 
of pixel signal to noise (apart from a few outliers 
due to bright sources) when the variance is calcu- 
lated in this manner. 

Pixel covariances were derived directly from the 
image data. We calculated the average covariance 
for all pixel pair configurations (up to a separa- 
tion such that the covariance was negligible) over 
image areas away from bright sources. Avoiding 
areas with bright sources helps minimise the over- 
estimation caused by the instrumental PSF. Nev- 
ertheless, as we discuss in Section[3j the covariance 
is still overestimated by a small amount, meaning 
that our quoted significances are conservative. 

3. Analysis 

Rather than follow the procedure of stack- 
ing small sections of the ima ge centred on th e 
source positions (see, for e.g., iDole et al. 2006), 
we opt for the method used in our earlier work 



(|Dye et al. 1120061 ) whereby flux is measured di- 



rectly from the image at each source position. The 
catalogue of sources is offset by varying amounts 
on a 2D regular grid and the flux summed over 
all sources at each offset. The result is an 'offset 
map' that gives an indication of how well aligned 
the sources are with respect to the image and the 
significance o f the stacked flux (see Figure [2|) . As 



we showed in lDve et al. (2006), if the sources are 



properly aligned with the image, then the correct 
stacked flux is that at the origin of the offset map, 
not necessarily at the peak which may be slightly 
offset from the origin. 

( 20061) 



The data stacked in iDve et al 



were 

SCUBA maps with each pixel value representing 
the total flux a point source would have if located 
within that pixel. The stacking therefore simply 
took the sum of all map pixel values at the source 
positions. In the current work, the MIPS data 
output by the pipeline adhere to the usual opti- 
cal convention whereby a pixel holds the flux re- 
ceived solely by that pixel. A source's total flux 
is therefore the sum of flux in all pixels belonging 
to the source. To convert the MIPS data into the 
convention used by the SCUBA data in prepara- 
tion for stacking, we convolved the images with 
a circular top-hat, then multiplied them by the 
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aperture correction corresponding to the top-hat 
radius. The convolution was carried out at the 
original pixel scale of each image, so to prevent 
aliasing effects, pixels around the top-hat circum- 
ference were weighted by their interior fractional 
area. 

Our choice of a circular top-hat instead of the 
more conventional instrument PSF was based on 
the fact that the MIPS PSF varies between sources 
and between images. We created simulated MIPS 
images of a point source, varying the asymmetry 
and size of the image PSF compared to the fiducial 
model PSF in each case. We found that the error 
in the total source flux measured by convolving 
with the fiducial PSF rises more quickly with in- 
creasing PSF asymmetry and size than measured 
by convolving with a circular top-hat having an 
aperture correction matched to the fiducial PSF. 
The MIPS data handbook recommends that the 
PSF should be determined directly from bright 
sources in the image, but since our image has no 
sufficiently bright sources, this was not possible. 

With this in mind, we chose top-hat radii of 
r = 18" and r = 40" for the 70 and 160/im 
images respectively. Instead of using the corre- 
sponding aperture corrections from the MIPS data 
handbook, we computed our own to ensure con- 
sistency with our top-hat convolution. For each 
wavelength, we took the in-orbit PSlH, binned it 
to the relevant image pixel scale, then computed 
its product with the edge-weighted top-hat to give 
the fraction of flux contained within the top-hat 
and hence the aperture correction. For the 70/im 
data, we measured an aperture correction of 1.63 
for the r = 18" top-hat and for the 160/im data 
with the r = 40" top-hat, an aperture correction of 
1.53. These are ~ 5% smaller than the low tem- 
perature aperture corrections given in the MIPS 
handbook. 

In this paper, we quote an average stacked flux 
per source, / = Y>f =1 fi/N and an average in- 
verse variance weighted flux per source, f w = 
T>i(fia~ 2 )/Y,ia~ 2 . Here, /j and Oi are respectively 
the flux and la uncertainty on the top-hat con- 
volved image pixel populated by source i. Com- 
parison of the average flux with the weighted aver- 
age flux gives an indication of whether the stacked 
signal is dominated by a minority of high signif- 



icance sources. Simple error propagation shows 
that ai is given by 



(i) 



where U is the value of the top-hat function in 
pixel i and Cjk is the covariance in the original, 
unconvolved image between pixels j and k. As 
explained in Section 12.21 the variances, i.e., diag- 
onal terms of Cjfc, come from the variance image, 
computed for each image pixel from the coverage 
map. However, the off-diagonal terms are the co- 
variances averaged over the whole image, so that 
any pixel pair jk with the same separation vector 
are assigned the same covariance. 

To verify our treatment of errors, for each of 
the 70 and 160/mi data, we fitted a Gaussian to 
the distribution of flux significance in the origi- 
nal unconvolved images, and to the distribution 
of the significance, /i/c,-, for the convolved im- 
ages. In the latter case, we first included, then 
omitted the off-diagonal elements in the covariance 
matrix. With the original 70 and 160/xm images, 
the Gaussian fit had unit standard deviation as 
expected. With the convolved images and the off- 
diagonal covariance terms included, the standard 
deviation for the 70/im image was approximately 
0.95 and for the 160/im image, 0.90. However, in- 
cluding only the variance terms gave a standard 
deviation of 1.79 for the 70/im data and 2.78 for 
the 160/im data. This test confirms two facts: 1) 
If covariance is not allowed for, the stacked flux 
error is underestimated by ~ 45% at 70/im and 
~ 65% at 160/xm. 2) Our measurement of covari- 
ance is slightly overestimated, presumably due to 
the MIPS PSF, giving rise to a conservative 5-10% 
underestimate of the stacked flux significance. 

Finally, source confusion due to the large 70 
and 160/im MIPS PSF must also be accounted for 
in the stacking. With the nodde d and chopped 
SCUBA data of iDve et al.l (|2006h . this could be 



2 provided at http://ssc.spitzer.caltech.edu/mips/psf.html 



neglected because the beam and hence the map in 
these data had an average of zero. With the MIPS 
data, this is not so. The average stacked flux per 
source must therefore be corrected by subtracting 
off the average of the convolved image then divid- 
ing the result by the factor (1 — B e /A), where B e 
is the effective area of the PSF of the convolved 
image and A is the image area (see Appendix |A"| . 
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4. Results 

To investigate the contribution from bright, di- 
rectly detectable sources in the MIPS images to 
the average stacked flux, we carried out two stacks 
per image, one leaving the image unaltered and a 
second with all > 3cr sources removed. At 70/im, 
there are six > 3a sources, whereas at 160/im, 
there are three (see Figure [T]). Sources were re- 
moved by subtracting the in-orbit PSFs from the 
images at the position of each source, scaled to 
match the integrated source brightness. 

4.1. Stacking the full SCUBA, 24/im and 
8/im catalogues 

The results of stacking all (i.e., not selected by 
redshift) SCUBA, 24/im and 8/im sources onto the 
MIPS data are given in Table[T] Offset maps show- 
ing the average weighted flux per source for each 
combination of MIPS image and source list are 
also plotted in Figure [2] Errors in both the ta- 
ble and the maps include the uncertainty of the 
calibration on the 70 and 160/im data. 

Apart from a single case, i.e., the instance 
in which SCUBA sources were stacked onto the 
160/im MIPS image with > 3a sources removed, 
every stacking combination results in a significant 
detection of far-IR flux. All peaks in the offset 
maps are well aligned with the origin. This in- 
dicates that all data are properly aligned, as we 
expected since all Spitzer data are tied to 2MASS 
(Two Micron All Sk y Survey; ICutri et al. 1 120031) 

(|2006l ) that the 



and we showed in iDve et al 



SCUBA data are well aligned with the Spitzer 
data. The more extended nature of the peaks in 
the 160/im offset maps is a reflection of the broader 
PSF at this wavelength (40"FWHM, compared to 
18"at 70 urn). 

The differences between the 3a source-subtracted 
and unmodified stacks show that at both 70 and 
160/im, the > 3a sources account for approxi- 
mately 50% of the average flux per source, across 
all three source populations. In every case, the av- 
erage flux per SCUBA source is highest, followed 
by the average flux per 24/im source then per 8/im 
source. This is not surprising; the 70 and 160/im 
data are sensitive to the same dusty population 
of sources as SCUBA, whereas the 8/im data are 
also sensitive to distant older stellar populations. 
Also, the fact that there are more objects in the 



24 and 8/im catalogues brings the average flux 
down because on average, these sources will sam- 
ple more image noise than areas of significant 70 
and 160/im emission. 

Dole et al. I (2006) stack 24/im sources with 



fluxes > 60/iJy onto MIPS data to measure a flux 
of 0.138 ±0.024 MJy/sr at 70/im and 0.571 ±0.123 
MJy/sr at 160/im. Using our MIPS data with 
all > 3cr sources removed, we find a lower con- 
tribution of 0.070±0.010 MJy/sr at 70/im and 
0.36±0.09 MJy/sr at 160/im. However, a con- 
tribution of 0.103±0.019 MJy/sr at 70/im and 
0.87±0.16 MJy/sr at 160/im is made by the 8/im 
sources, again having removed the > 3a sources. 
Within the errors and including the fact that our 
data are more prone to cosmic variance (see be- 
low) being ~ 80 times smaller in are al coverage. 



our re sults are consistent with those of lDole et al. 
(|2006h . 



To estimate of the effects of cosmic variance on 
our results, we divided the data into two approxi- 
mately equal areas and then repeated the stacking 
with the halved data. This was performed twice, 
firstly splitting by the median source RA of each 
source catalogue, then by the median Declination. 
The la variation in the spread of the resulting 
stacked 70/im flux was found to be ~ 30% and 
the variation in the 160/im flux, ~ 20%. Since 
this is an estimate of the variance between fields 
half the size, the variance between fields of the full 
size, discounting clustering effects, will be a factor 
of y/2 smaller, i.e., - 20% at 70/im and - 15% 
at 160/im. As this serves merely as an order-of- 
magnitude estimate of the cosmic variance, we do 
not include it in any of the errors quoted in this 
paper. 

Contribution of sources to the CIB 

The total contribution of the three different 
source populations to t he CIB is given i n Ta - 



ble [2] By extrapolation, IChapman et al. ] (|2005l) 



estimated an upper limit on the contribution of 
> lmJy SCUBA sources detected at 850/im to the 
CIB emission at 160/im of <0.04 MJy/sr. Despite 
our CUDSS sources having a brighter sensitivity 
level of 3.5mJy, we measure a higher contribu- 
tion at 160/im of 0.125±0.040 MJy/sr, without re- 
moving any bright MIPS sources, or 0.060±0.042 
MJy/sr having removed all > 3a sources. Al- 
though these measurements have large uncertain- 
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Fig. 2. — Offset maps of average weighted stacked flux per source (f w ) for all combinations of MIPS image and 
sources stacked. Columns from left to right correspond tp SCUBA sources, 24/im sources and 8/im sources. First 
and second rows correspond to 70/im MIPS data including then excluding > 3a sources respectively. Similarly, third 
and fourth rows correspond to 160/im MIPS data including then excluding the > 3<j sources. Contours start at 2a 
significance and increase by 0.5<r intervals. Average fluxes for each case are given in Table[T] Significances account for 
the MIPS 70 and 160/im calibration uncertainty. The number of sources stacked for each of the different combinations 
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MIPS Data 


850Mm SCUBA Sources 


24/im Sources 


8/im Sources 


70/im 

70/im - 3ct sources 
160/im 

160/im - 3cr sources 


3.63±0.77 (3.45±0.80) 
2.51±0.77 (2.58±0.78) 
19.9±6.4 (16.9±6.6) 
9.5±6.6 (8.1±6.7) 


2.04±0.25 (2.10±0.25) 
1.64±0.24 (1.65±0.24) 
15.0±2.4 (14.9±2.4) 
8.8±2.2 (8.6±2.2) 


0.95±0.12 (0.83±0.12) 
0.64±0.12 (0.57±0.12) 
8.5±1.1 (8.2±1.1) 
5.6±1.0 (5.4±1.0) 



Table 1: Average weighted flux per stacked source (f w ) in mJy for the MIPS 70 and 160/im images, including 
and having subtracted > 3a sources. These correspond to the flux at (0, 0) in the offset maps shown in Figure 
[2j Quantities in parentheses are the average unweighted fluxes /. All errors include the MIPS 70 and 160/rni 
calibration uncertainty. 



ties, they suggest the possibility of a somewhat 
larger SCUBA source contribution to the 200/im 
CIB peak than previously thought. 

Table [2] shows that at both 70 and 160/im, 
the SCUBA sources make the lowest total contri- 
bution, followed by the 24/xm sources and then 
the 8/xm sources with the highest contribution. 
This is an important result; sources on the shorter 
wavelength side of the peak in the CIB resolve 
more of the CIB at 70 and 160/im, and there- 
fore most likely at the 200/im peak itself, than 
the SCUBA sources on the longer wavelength side. 
This is almost entirely due to the differing sensi- 
tivities of the source populations used for stack- 
ing. In terms of the efficiency of resolving the 
bulk of the CIB emission, the 8 and 24/im source 
population are more favourable than the SCUBA 
sources. This is not surprising because SCUBA 
surveys typically find 0.4 sources per arcmin 2 for 
each hour of observation whereas Spitzer surveys 
find ~ 130 8/im sources per arcmin 2 for each hour 
of observation. Of course, in the context of the 
present study, SCUBA's time would be more effi- 
ciently used by computing the cross correlation of 
the 850/im maps with the Spitzer images, rather 
than merely stacking at the positions of signifi- 
cant SCUBA sources. We will carry out this cross 
correlation in future work. 

Expressing the absolute contributions in Table 
Has a fraction of the CIB is somewhat difficult 
due to the uncertainty in the background flux at 
70 and 160/im (primarily because of differing esti- 
mates of foreground contamination). In fact, there 
are no direct measurements at these specific wave- 
lengths. The most reliable measurements close 
to 160/im are those at 140/iin made by the Dif- 
fuse Infrared Background Experiment (DIRBE) 
on board the Cosmic Background Explorer. De- 



pending on the calibration used, the DIRBE re- 
sults give a flux of eith er 1.17 ± 0.32 MJy/ sr or 



0.70 MJy/sr at 140 /im dHauser et al. I ll998h. al- 
though as noted bv lDole et al. H2006), th e zodi 



acal cloud colours of iKelsall et al. I ( 19981 ) imply 
that a further 0.14 MJy/sr should be subtracted 
from these numbers. Taking the average of both 
calibrations and subtracting 0.14 MJy/sr gives a 
flux of 0.80 MJy/sr. This can be extrapolated 
to give an approximation of the flux at 160/im 
of 0.99 MJy/ s r usin g the SED fit to the CIB by 



Fixsen et al. 



( 1998h . DIRBE also provided esti- 
mates of the CIB at 60/im which are a useful con- 
straint on our mea surement of the background at 



70/im with MIPS. iFinkbeiner. Davis fc Schlegel 



(2000) placed an upper limit on the CIB at 60/im 
of 0.56 ± 0.14 MJy/sr. Mor e recently, this limit 
was re duced to 0.3 MJy/sr bv lDwek fc Krennrich 
(I2OO5I) . 



In terms of lower limits on the CIB at these 
wavelengths, stacking analyses are currently the 
most stringent. By e ffectively extrapola ting their 
24/im number counts, Dole et al. (|2006f ) currently 
provide the highest lower limits on the CIB of 



MIPS Data 


SCUBA 


24/im 


8/tm 


70/im 
70/im - 3cr 
160/im 
160/im - 3cr 


0.021±0.005 
0.014±0.004 
0.125±0.040 
0.060±0.042 


0.087±0.011 
0.070±0.010 
0.62±0.10 
0.36±0.09 


0.152±0.019 
0.103±0.019 
1.31±0.17 
0.87±0.16 



Table 2: Contribution of sources to the CIB at 70 
and 160/im in units of MJy/sr, computed from the 
average weighted flux per source given in Table 
[TJ The coverage for each source catalogue is as 
follows: 41 arcmin 2 for 850/im SCUBA sources, 
49 arcmin 2 for the 24/im sources and 59 arcmin 2 
for the 8/im sources. 



S 



0.17±0.03 MJy/sr at 70/im and 0.71±0.09 MJy/sr 
at 160/im. Using these and taking the DIRBE 
measurements quoted above as upper limits, we 
can estimate the range in fractional contribution 
that our strongest contributors, the 8/im sources, 
make to the CIB. The upper limit to t he CI B at 
60/im imposed bv lDwek fc Krennrich ( 2005) and 



the extrapolated lower limit of iDole et al. 



at 70/im indicates that our 8/im sources resolve 
~ 35 — 75% of the background at 70/mr. Simi- 
larly, taking the DIRBE extrapolation to 160/im 
as an upper lim i t and the extrapolated lower limit 
of IDole et all (|2006h at 160/im implies that ~ 
90 — 100% of the 160/im background is resolved 
by the 8/im sources. 

To what extent do the additional 8/im sources 
not detected in the 24/im data emit at far-IR wave- 
lengths? This can be very crudely estimated by 
calculating the number of 8/im sources that would 
be required to give the same measured CIB con- 
tribution but assuming each source has a constant 
flux equal to the corresponding average 24/im 
source flux. Here, the assumption is made that 
all the 24/im sources (95% of which are detected 
at 8/iin) contribute to the far-IR flux. This simple 
calculation, shows that ~ 20% and ~ 50% of the 
additional 8/im sources at 70 and 160/im respec- 
tively would have to contribute in that case. This 
is a conservative estimate because in reality, the 
average far-IR flux of the additional 8/im sources 
will be lower than the average flux of the 24/im 
sources due to the increased sensitivity of IRAC 
at 8/im. 

4-1.2. Average 8 & 24^m source SEDs 

InlDve" et al. ] (|2006h we measured the average 
450 and 850/iin flux per 24 and 8/im source. Com- 
bining these measurements with the average 70 
and 160/im flux per 24 and 8/im source deter- 
mined in the present work gives four data points 
each to which we can fit average SEDs. Figure [3] 
shows the results of fitting the grey-body function 
Ai/PB^, T) to these average fluxes, where A is a 
normalisation constant and B is the Planck func- 
tion. In the fit, the parameters A, f3 and T were 
allowed to vary and we redshifted the function to 
the median redshift of our sample, z = 1.0. 

For the 24/im sources, the best fit is achieved 



the 8/im sources with j3 = 1.54j^o|g and T = 
38.8l^K (la errors quoted). The dependence 
of these fitted parameters on the median redshift 
is such that a change in redshift Az produces a 
change in T given by AT ~ 19Az for both 24 
and 8/im sources, while /3 has absolutely no depen- 
dence. The temperature of our average sources is 
consistent with temperatures of submm galaxies 
found in the local universe, e.g., iDunne fc Eales 



(|200lf l who measure T=(36± 5)K. Whether there 
is consistency with submm sources in the high 
redshift Universe is less clear. The sample of 73 
SCUBA sources with a median redshift of 2.3 of 
Chapman et al. I {2005) has a median temperature 
of T mc( j ~ (36 ± 7)K, consistent with our values. 
However, the sample of 1 SCUBA sources w ith a 
median redshift of 1.7 of Pope et al. (j2006f ) has 



a lower median temperature of T me( j ~ 30K. If 
a discrepancy exists, then it could be explained, 
at least in part, b y the selection effect noted 
by Chapman et al. I (2005); surveys like that of 
Pope et al. ( 20061 ) requiring a submm and radio 
detection are biased toward colder sources. 

The normalisation of both SEDs confirms 
that the average 24 and 8/im source detected 




o.oi 



100 



1000 



Observed wavelength (jim) 



with (3 = 2.05±o ^ and T 



39.7±^K and for 



Fig. 3. — SED fits to the average 8/im (black line, 
open circles) and 24/im (grey line, filled circles) source 
fluxes for the median source redshift of z = 1. Data 
points at 70 and 160/im are taken from the present 
study and those at 450 and 850/im from iDve et al. I 
( 2006). Both SEDs show that the average 8 and 24/im 
sources detected by Spitzer in this dataset are border- 
line ULIRGs. 
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by Spitzer in the current sample is a border- 
line ultra- luminous infrared galaxy (ULIRG). 
To demonstrate this, we use th e defini tion of 
IClements. Saunders fc McMahonl (| 19991 ) that 
stipulates a ULIRG must have a luminosity of 
at least 2.5 x 1O 11 L measured at 60/mi by the 
Infra-Red Astronomical Satellite (IRAS). The 
rest-frame 60/iin flux computed from our best fit 
SEDs is 2.3 x 1O 11 L0 for the average 24/im source 
and 1.2 x 10 11 L n for the averag e 8/xm source, in 



good agreement with lDve et al. ] (|2006l) 



An interesting question is how do our aver- 
age 24 and 8/im sources compare to the aver- 
age SCUBA sourc e detected by other studies? 



Pope et al. 



( 2006h define the quantity L m as 
the integral of flux in the wavelength range 8 - 
1000/im. Their sample of 10 SCUBA sources has 
a median Lir of 6.0 x 10 12 L Q . Similarly, the me- 
dian valu e of Ltr for the s ample of 73 SCUBA 
sources of IChapman et al. (2005) is 8 x 10 L Q . 
In comparison, our average 24/im SED gives Lm, = 
5.8 x 1O 11 L and our average 8/im SED Ljr = 
3.5 x 10 11 Lq. The average 24 and 8/im source in 
our sample is therefore ~ 10 times fainter than the 
average SCUBA source detected by the previous 
two studies. 

4.2. Stacking the 24 and 8/im sources by 
redshift 

In this section, we consider the contribution 
from the 24 and 8/im sources to the CIB at 70 and 
160/mi as a function of redshift. Using the photo 



metric redshifts already determined in lDve et al. 



( 20061 ). we divided the sources equally into redshift 
bins of varying width. The source redshifts extend 
up to z ~ 4 (see Figure 3 of Dye et al. 1120061 ). Bins 
were chosen to be large compared to the average 
redshift uncertainty but small enough to give rea- 
sonable resolution, hence the 24/im sources were 
divided into 5 redshift bins and the 8/im sources 
into 6. Approximately 10% of sources with un- 
determined redshifts (due to their sparse optical 
photometry) were omitted from the analysis of 
this section. This therefore gives ~ 30 objects per 
24/im bin and ~ 120 objects per 8/im bin. 

Figure|3]shows how the weighted average 70 and 
160/im flux per source varies with redshift. The 
plots show this variation both having removed the 
> 3(7 sources in the MIPS images and with them 
left in place. At 70/iin, the effect of removing the 



3tr sources has less effect than at 160/im. Also, at 
70/im, the flux is dominated by 8/im sources ly- 
ing at lower redshifts. Dividing the 8/^m sources 
into two populations segregated by the median 
redshift, z = 1.0, the low redshift population ac- 
counts for (79 ± 10)% of the total 70/im emission 
(having removed the > 3er sources) from the 8/im 
sources. In comparison, the 70/im emission per 




Fig. 4. — Variation of average weighted 70/im flux 
(top) and 160/im flux (bottom) of Spitzer 24 and 
8/im objects binned by redshift. The median red- 
shift in each b in is plotted. Re dshifts extend up 
to z ~ 4 (see iDve et al. 2006 ) and are divided 
equally between bins. In both plots, the continu- 
ous grey and black lines correspond to the 24 and 
8/im sources respectively, stacked onto the MIPS 
images without any sources removed. The dashed 
lines show the average weighted flux having re- 
moved the > 3cr sources from the MIPS images. 
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24/im source is more evenly spread in redshift, the 
low redshift population accounting for (51 ± 10)%. 
At 160/im, the low redshift 8/im sources contribute 
(52 ± 14)% of their total and the 24/im sources 
contribute (42 ± 15)%, having removed all > 3a 
160/im sources. 

The differences between the 70 and 160/im plots 
in Figure 2] are very well explained by a single av- 
erage source SED consistent with the fitted aver- 
age SEDs derived in section l4\1.2| To demonstrate 



this, we took an SED from Dale & Helou 



( 20021) 

corresponding to a dust temperature of T=40K to 
match our average SEDs. Since this SED extends 
into the optical and models typical mid-IR spec- 
tral features due to dust, a realistic prediction of 
the 70 and 160/im flux of a source given its red- 
shift and 8 or 24/im flux can be made. In this 
way, using our 8 and 24/im source catalogues, we 
computed a prediction of the variation of 70 and 
160/im flux with redshift. 

The results of this analysis are shown in Fig- 
ure [5l There is good agreement between our 
measured variation and the predicted variation. 
We reproduce the total flux (summed over all 
sources) within the errors and also the observed 
trends. Most notably, we reproduce the decline in 
the 70/im emission from the 8/tm sources within 
< z < 1 as the peak of the SED is redshifted 
out of the 70/im band. This explains why the ma- 
jority of 70/im emission is observed from the 8/im 
sources lying at z < 1. The prediction degrades 
quickly if a cooler or warmer SED is used; with 
a 35K or 45K SED, the total predicted flux is in- 
consistent with the total measured. The fact that 
a single SED can be used to fit the observed flux 
over such a wide range of redshifts implies that 
only a small amount of source evolution must have 
occured during that time. 

To assess the effects of cosmic variance on the 
results of this section, we repeated the previous 
exercise of dividing the data into halves and re- 
stacking. We found that the major trends are ro- 
bust, i.e., the decline in 70/im flux from the 8/im 
sources over the interval < z < 1 and that the 
other combinations remain consistent with little 
or no variation with redshift. However, the large 
spike in 160/im flux seen from the 24/im sources 
at z ~ 1.5 (without having removed the > 3a 
sources) is not robust and therefore presumably 
an effect of cosmic variance. 



5. Summary and Discussion 

In this paper, we have quantified the contri- 
bution of flux to the CIB at 70 and 160/im from 
SCUBA sources and 24 and 8/im Spitzer sources in 
the CUDSS 14hour field. By stacking flux at the 
position of the different sources, we have found 
that the SCUBA sources make the highest contri- 
bution per source and that the 8/tm sources make 
the lowest. Conversely, the opposite is true of 
the total contribution from all sources, leading to 
the conclusion that the bulk of the CIB is most 



70|im 




0.5 1 1.5 



2.5 




Fig. 5. — Comparison of the measured (dashed 
lines) and predicted (continuous lines) MIPS 70 
and 160/im flux from the 8/im (black) and 24/im 
(grey) sources using a T=40K SED taken from 



Dale fc Helou I (|2002f ) . The measured data are the 
stacked fluxes in Figure 2] with > 3a sources re- 
moved. 
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efficiently resolved by sources detected at wave- 
lengths shorter than the peak of the CIB emission 
at ~ 200/im. 

Our stacking suggests a somewhat larger con- 
tribution from the CUDSS SCUBA sources to 
the 200/im CIB peak than previously thought. 
Chapman et al. I (|2005f ) estimated an upper limit 
on the contribution of > lmJy SCUBA sources 
detected at 850/im to the CIB emission at IGO/mi 
of <0.04 MJy/sr. Despite our SCUBA sources 
having a brighter sensitivity level of 3.5mJy, we 
measure a contribution at 160/im of 0.125±0.040 
MJy/sr, without removing any bright MIPS 
sources, or 0.060±0.042 MJy/sr having removed 
all > 3cr sources. 

Since measurements of the CIB at 70 and 
160/im are presently very uncertain, the fractional 
contribution to the CIB made by our sources can 
only be expressed within a range set by present 
upper and lower limits. Using the DIRBE esti- 
ma tes as upper limits and the lower limits set 
by iDole et al. (2006), our strongest contribu- 
tors, the 8/im sources, resolve somewhere be- 
tween ~ 35 — 75% of the background at 70/im 
and - 90 - 100% at 160/im. 

By combining our results in the present work 
with our pre vious stacking of 450 and 850/im 
SCUBA flux (|Dve et al. I l2006h . we have estab- 
lished that the 8 and 24/im sources detected by 
Spitzer are on average borderline ULIRGs. The 
average source we detect is ~ 10 times fainter 
than the a v erage SCU BA source detec ted by 



tnan tne a v erage am riA source detec t ed Dy 
Pope et al. I (|2006h and IChapman et al. I (2005) 



integrating flux over the wavelength range 8 - 
1000/im. Furthermore, the temperature of ~ 40K 
of our average source is consistent with tem- 
peratures of su bmm galaxies found i n the local 
Universe (e.g., Dunne fe Eales I l200lh and with 
the median temperature of (36 ± 7)K of SCUBA 
sources in the h i gh re dshift Universe found by 



Chapman et al. I (|2005l ). However, our average 
source is warmer than the median temperature of 
T ms d ~ 30K of 10 SCUBA sources measured by 



Pope et al. 



(|2006h . 

Using photometric redshifts assigned to the 8 
and 24/xm sources, we have investigated how the 
contribution of 70 and 160/im flux to the CIB 
varies with redshift. We have found that the 8/im 
sources at low redshifts, z < 1 (accounting for 
half of them), are the strongest contributors to 



the CIB at 70/zm, their flux amounting to ~ 4 
times that of the z > 1 sources. The 70/im emis- 
sion per 24/im source as well as the 160/im emis- 
sion per 8 and 24/im source is consistent with an 
even distribution over redshift. This verifies the 



result of IDole et al. I (2006) that the majority of 
the emission at 70 and 160/im from 24/im sources 
must come from a redshift of z ~ 1 where the 
redshift distribution of these sources peaks. We 
have shown how this distribution can be repro- 
duced from our observed 8 and 24/im catalogue 
of fluxes and redshifts using a single non-evolving 
source SED with a dust temperature of 40K. 

As a concluding remar k, this study and sim- 
ilar recent studies (e. g., Serjeant et al. 20041 



Knudsen et "aLl2005t[Dve et al. ll2006tlWang. Cowie fc Barger 



2006) indicate that the CIB is not predomi- 
nantly due to a rare population of exception- 
ally luminous submillimeter sources as hinted 
at by early SCUBA observations, but that a 
much more numerous galaxy population of mod- 
est average luminosity is responsible instead. 
However, we have shown here that the SCUBA 
galaxies probably do make a larger contribu- 
tion than previously thought, although much 
larger numbers of SCUBA sources such as those 
of the SCU BA half degree extr agalactic survey 



(SHADES; iMortier et all 12005 ) or those result- 
ing from future SCUBA2 surveyCjwill be required 
to improve the precision of these measurements. 



3 see http://www.roe.ac.uk/ukatc/projects/scubatwo 
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A. Flux boosting correction 

In the following, it is assumed that the MIPS image has been prepared such that the value of any one pixel 
gives the total flux a point source would have if located within that pixel. Since this preparation inevitably 
involves convolution of the raw MIPS image with some kind of kernel, the profile of a point source will be 
the convolution of this kernel with the original image PSF. This resulting profile is referred to as the image 
'beam' hereafter. 

Suppose that T is the total number of sources being stacked onto the MIPS image and that a subset N of 
these are 'genuine' sources. A source is defined as 'genuine' if it has associated MIPS emission that makes a 
non-negligible contribution to the stacked flux. (In practice, one would derive a threshold flux that depends 
on the number of sources). The actual average flux per source is then simply the summed flux from the 
genuine sources divided by the total number of sources, 



1 N 

/actual = 7^ E fi > (Al) 



where fi is the flux in the MIPS image pixel at the position of the genuine source i. 

However, this quantity is overestimated if one naively adds the flux in the MIPS image at the positions of all 
T sources for two reasons. Firstly, the T — N sources without associated MIPS emission (the 'contaminating' 
sources) sample flux from the genuine sources since some will happen to lie within genuine source beams. 
The extra flux sampled on average from a genuine source with flux /j by the contaminating sources is 

ft = 2™ c /, J rb(r)dr = n c B e f i} (A2) 

where b(r) is the radial beam profile scaled to have a peak height of unity, n c is the number density of 
contaminating sources and B e defines the effective beam area. Secondly, the genuine sources themselves 
sample emission from neighbouring genuine sources when their beams overlap. Similar to the contaminating 
sources, the extra flux sampled on average from a genuine source with flux /, by its neighbouring genuine 
sources is 

// = n' g B e f z , (A3) 

where n' g is the number density of the neighbouring N — 1 genuine sources. The average MIPS flux per 
source that is measured by summing the flux at all T source positions is therefore 

N 1 l / T3 N 

1 + n,iL 



measured 



i i 

= (l + n;B e )/ actual , (A4) 

where n' t = n' g + n c is the number density of T — 1 sources. The measured average flux is therefore the actual 
average flux boosted by a factor of (1 + n' t B e ). 

If the MIPS data are properly normalised, (so that there is no net positive emission from artifacts, 
systematic effects, etc.) then the average of the image over its area A is 

- B N n t B N - 

/image = E & ~ J 1 E ^ = n *^ e /actual ' (A5) 

i i 

having used the fact that the number density of all T sources is n t — T/A. Subtracting equation (|A5[) from 
equation (|A4j) and rearranging gives 

/actual = (/measured — /image) 0- ~ Be/ A) , (A6) 
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hence the actual average flux per source can be obtained by subtracting the image average and dividing by 
the factor (1 — B e /A). 

In the above, it has been assumed that source positions are a random sampling of a uniform distribution. 
In reality, the sources will exhibit a degree of clustering. Clustering of the genuine sources causes a positive 
bias of the average flux compared to the non-clustered assumption, whereas clustering of the contaminating 
sources on average has no effect. 
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